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Abstract  of  Technical  Progress 

In  earlier  work  the  necessary  conditions  of  optimality  were  derived  for  a 
problem  of  minimum  miss-distance  guidance  of  air-to-air  missiles.  The  model  was 
based  upon  nonlinear  translational  equations  of  motion.  The  solution  of  the 
necessary  conditions  requires  a  solution  of  a  two-point  boundary-condition  problem. 
Two  methods  proposed  for  the  latter  solution,  an  elliptic  integral  method  and  a 
series  technique,  were  studied  and  both  methods  were  rejected  in  favor  of  a  pro¬ 
cedure  based  upon  the  quasilinearization  method.  The  latter  requires  fewer  assump¬ 
tions  and  exhibits  excellent  convergence  prooerties. 

In  order  to  remove  the  numerical  integration  problem  ana  to  simplify  the  linear 
two-pomc  boundary— cond i t ion  problem  associated  with  quasiiinearizacion,  the  regular 
method  was  modified,  three  alternative  techniques  being  derived,  and  a  technical 
report  was  written  which  discusses  the  convergence  properties  and  accuracy  of  the 
three  modified  quasiiir.earization  methods  applied  to  two-point  boundary-condition 
problems  m  general.  . 

One  of  the  latter  mecnoas  was  applied  to  the  missile  guidance  problem  and  com¬ 
pared  to  a  linear  methoa  of  guidance.  The  comparison  was  favorable  in  the  Cuies  of 
missile  encounters  considered.  A  second  technical  report  was  prepared  which  des¬ 
cribes  the  application  and  the  numerical  results. 
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Technical  Report 


The  principal  objective  was  the  development  of  an  optimal  guidance  method 
fo;  air-to-air  missiles.  The  proposed  model  was  based  upon  nonlinear  cranslational 
equations  of  motion  and  a  quadratic  performance  index.  Availability  of  a  predicted 
time-history  of  the  target's  posicion  was  assumed. 

In  an  earlier  investigation  the  necessary  conditions  of  optimality  were  derived 
and  two  methods  were  proposed  for  solving  these  conditions.  They  were:  An  elliptic 
integral  method  and  a  series  method.  These  were  to  be  placed  in  a  form  suitable  for 
computation  and  applied  to  the  solution  of  the  necessary  conditions,  which  reduce  to 
a  two-poir.c  boundary-condition  (TP3C)  problem  associated  with  a  system  of  nonlinear 
ordinary  differential  equations.  As  far  as  possible  the  following  extensions  ;among 
others)  to  the  elliptic,  integral  and  series  methods  were  to  be  made: 

(a)  Introduction  of  variable  missile  velocity  magnitude; 

(b)  Extension  into  three  dimensions; 

<z)  Allowance  for  variations  of  more  than  90"  in  the  flight  direction  angles; 

(d )  Inclusion  of  the  case  of  control  angles  which  do  nor  vary  monotonicallv 
with  time; 

(e)  Inclusion  of  an  improved  method  for  calculating  Cime-to-go; 

(f)  Placing  bounds  upon  the  control  variables. 

After  a  study  of  the  methods  based  upon  elliptic  integrals  and  series,  it  was 
decided  that  a  third  method  based  upon  quasilinearization  (or  the  generalized  Nr  ’ton- 
Raphson  method)  would  be  a  more  powerful  computational  device  (as  far  as  numerical 
convergence  to  a  solution  to  the  TPBC  problem  is  concerned)  and  a  more  flexible  tool 
in  regard  to  the  implementation  of  extensions  (a)  -  (£)  given  above.  On  consulting 
with  the  Program  Manager,  it  was  decided  that  the  quasi  linear izat ion  method  would  be 


used  rather  than  the  methods  originally  proposed  for  study. 
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However,  the  usual  quasilinear izacion  algorithm  Is  involved,  computationally 
speaking,  because  it  requires  many  numerical  integrations  of  systems  of  differential 
equations  and  the  solution  of  many  linear  TPBC  problems.  With  this  in  mind  the 
following  enclosed  technical  report  was  prepared: 

[1]  "Modified  Methods  of  Quasilinearization". 

The  latter  report  develops  modified  methods  of  quasilinearizacion  .n  which  numerical 
integration  is  unnecessary  and  in  which  the  linear  TPBC  problems  to  be  solved  are 
reduced  in  number.  The  modified  quasilinearization  methods  are  applicable  to  any 
TPBC  problem.  A  shortened  version  of  the  report  has  been  submitted  for  possible 
publication  in  the  SIAM  Journal  of  N'umer leal  Analysis. 

One  of  the  modified  methods  of  quasilinearization  was  applied  to  the  air-to-air 
missile  problem,  the  details  being  given  in  the  second  enclosed  technical  report: 

[2]  "A  Guidance  Method  for  Air-to-Air 
Missiles  Using  Quasilinear ization" 

The  latter  report  contains  a  development  of  all  pertinent  equations  in  the  case 
of  encounters  in  three-d imens tonal  space.  It  shews  new  the  control  variables  may 
be  oounced  anc  gives  an  argument  (in  the  case  of  encounters  ir.  two-dimensional  space 
for  monotomc  variation  of  the  control  variable.  The  report  then  describes 
numerical  simulations  of  several  two-dimensional  missile-target  encounters.  The 
modified  method  of  quasilinearizacion  is  compared  to  a  linear  guidance  method,  the 
former  method  providing  more  accurate  commands  than  the  linear  method.  The  new 
method  allows  for  variable  nissle  velocity  magnitude,  determines  time-to-go,  and 
permits  large  changes  in  the  flight  direction  angles.  However,  the  largest  change 
which  has  been  simulated  is  90".  The  method  converged  to  the  optimal  solution  to 
the  nonlinear  guidance  problem  in  all  encounters  considered. 

The  only  technical  investigator  involved  in  this  study  has  been  Professor  J.  F. 
Andrus,  Ph.D.  in  Mathematics,  University  of  Florida,  June  1958.  His  thesis  was 
entitled  "Partially  Ordered,  Ideal  Preserving  Groups".  He  is  also  author  of  a 
number  of  papers  in  numerical  analysis,  optimal  control,  and  other  areas  of  applied 


matheoa  t ics . 


Modified  Quas ’ linearizat ion  Methods 


by 

J.  F.  Andrus 

Department  of  Mathematics 
University  of  New  Orleans 

I.  Introduction 

The  quasilinearization  (or  generalized  Newton  -  Raphson)  method  [2] 
has  proven  co  be  a  highly  effective  iterative  method  for  the  numerical  solution 
of  manv  two-point  boundarv-conditicn  (TPBC)  problems.  Frequently  only  a  rough 
approximation  is  required  to  start  the  convergence,  and  -  when  convergence  takes 
place  -  it  is  quadratic  m  character. 

'"-.ere  are  also  several  difficulties  associated  with  the  method.  They  are: 

(a;  The  differentia^  equations  must  be  linearized.  This  calls  for  iiff erer.tiatior. 

of  the  ng.hthand  :i__s  if  che  differential  equations. 

'b ;  After  each  iteration  the  approximate  solution  must  be  saved  for  use  in  the 

next  iteration.  This  presents  an  awkward  interpolation  problem.  Alternative¬ 
ly,  each  solution  mav  be  regenerated  by  means  of  integration  during  ail 
succeeding  iterations. 

.  During  each  iteration  ir.e  linear  differential  equations  must  oe  integrated 
''usually  numerically;  several  times  in  order  to  obtain  the  solution  to  a 
linear  two-point  boundary-value  (TP3V)  problem.  Sometimes  finite  difference 
techniques  are  used  instead  to  solve  the  TPBV  problem.  The  solution  is 
obtained  bv  solving  a  large  system  cf  linear  algebraic  equations. 

The  present  paper  contains  the  development  of  modified  quas 1  linear izac ion 
methods  in  vnich  the  problems  in  items  (b)and  (c;  nave  beer,  considerably  reduced. 
This  is  achieved  by  dividing  the  solutions  into  subarcs  and,  after  each  iteration, 
using  the  new  approximate  solution  to  obtain  constant  or  linear  approximations 
over  each  sub  interval .  The  equations  are  linearized  about  the  constant  or  linear 
approximations,  thereby  alleviating  the  problems  of  storage  and  interpolation. 

Moreover,  one  can  frequently  solve  the  linear  differential  equations  in 
terms  of  definite  integrals  or  even  in  closed  form  over  each  subinterval.  In 
such  cases  the  problem  of  numerical  integration  is  considerably  reduced. 
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Since  the  final  values  of  the  dependent  variables  on  each  subinterval  can 
be  expressed  linearly  in  terms  of  che  unknown  initial  values,  the  solutions  to 
the  complete  linear  TPSV  problem  can  be  easily  obtained  because  the  final  values 
of  the  dependent  variables  can  be  expressed  linearly  in  terms 

of  the  unknown  initial  values,  which  can  then  determined  bv  solving  a  system 
(usuailv  small)  of  linear  algebraic  equations.  (See  [1]  for  an  example.) 

It  is  argued  chat-under  some  reasonable  assumptions-  che  modified  methods 
of  quasilin.earization  will  converge.  Moreover,  convergence  becomes  quadratic 
in  nature  as  the  maximum  subinterval  size  is  decreased.  In  order  to  prove 
convergence  it  is  necessary  that  the  maximum  interval  size  be  chosen  sufficient! 
small.  It  is  also  proven  that  the  converged  approximate  solutions  to  the  -rigin. 
differential  equations  are  of  near  Iv  third— order  lccuracv  from  the  standpoint 
::  numerical  integration ;  in  the  case  or  the  constant  rerreser.ta  t  i  -ms  ■:  t.ne 
solutions  ever  subarts  anc  of  f ourth- order  iccuracv  m.  tn.e  cases  :f  toe  linear 
methods . 


Ir.e  would  ask  now 
third-order  accuracy. 

:  - 1  a  men  iron,  the  k-rn. 
i  certain  c  cn.s  t  ant  f  :r 
tier.  1 1  er  at  i  in  itself 


the  constant  solutions,  fer  eyamole,  cou.d  lead  to  nearly 
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:se  in  tn.e  k-i.  -tn.  iteraticn.  Therefore,  t.ne  quasi  line* 
on.tricutes  tc  th.e  .iccuracv  of  tn.e  tuner:  col  integral :  or.. 


II.  The  .'ms  i  linear  izac  ior.  Method 


ns.  ‘er  tc.e  T?EC  problem  consisting  of  the  inferential  equations 
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7h«>r»  is  Tt-  lssurir.-:  it  this  pnint  that  the  ao«i :  i  s  «*d  methods  of  tuasilmeu 

cion  will  corverge.  An  example  problem,  to  which  Method  1  has  been  applied,  is 
presented  in  Reference  [1]  in  full  detail.  In  the  latter  example  convergence  is 
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demonstrated.  N  the  convergence  problem  util  be  discussed  from  a  theoretical 
point  of  vi*>-  The  approach  will  be  to  transform  modified  quasilinearization 
problems  »«,  standard  problems  which  are  identical  in  every  respect  to  the 
modified  problems  except  for  certain  terms  which  approach  zero  as  t.ne  itera¬ 
tion  proceeds. 


The  problems  will  be  converted  to  standard  T?BC  problems  bv  referencing 
sub  intervals  to  a  single  interval  To  do  so,  let  it;  *  t;  -  t.  ,  and 
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In  the  case  of  Method  1  the  following  boundary  conditions  will  be  imposed: 
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As  N  »  •  and  max  ic^.  -  0,  the  srlutions  to  the  TPBC  problem  consisting 
of  equations  (3)  and  (4)  will  approach  the  solution  to  the  TTBC  problem  associatea 
with  aquations  1,2;  because  %  -n>  ^l°Lt)  . 

What  is  also  important  is  that  the  linearized  torn  of  equations  (3)  (including 
boundarv  conaitions)  will  be  the  same  as  that  of  equations  (2)  using  the  modified 
methods  of  quasiiinearioation  except  for  some  terms  unicn— under  certain  assump¬ 
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Also  to  be  considered  are  the  following  O.D.E.’s  oocained  by  setting 
u^  »  0  in  equations  (5): 
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equations  'a)  and  (5)*  are  equivalent  to  the  linearized  TPBC  probleits  associated 

»ith  the  three  modified  cuasiiinear iZ3tion  methods.  (It  is  assumed  that  all  end 
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were  to  be  izr.ored,  equations  (2)  anc  (1)  would  lead  to  the  same  quas ilineariza- 
ticr.  algorithms  as  would  the  nodi:  led  methods  applied  to  equations  1}  and  the 


conditions,  as  well  as  the  O.D.E.’s  nave  been  linearized.) 


original  boundary  conditions. 


„i  wii.  oe  assumea 


nat  it  :cr 


the  initial  approximation 


k*  ( •  \  r->t  i)  a  (i), 

i  *  )  is  sufficiently  close  to  the  true  solution  ,  z_  )  to  tne 


V  V 
'  3 

TPBC  problem  consisting  of  equations  (h)  and  \i),  then  the  standard  quasilir.eariza- 
:izn  method  a.  plied  to  the  latter  problem  will  converge  quadrat icallv  tc  the 
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differencial  equations  (5;*  and  the  boundarv  conditions  (4). 
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In  the  cases  of  Methods  2  and  3.  3  and  for  i  *  1,2,...,N  are  also  to  be 

-  .  .  uuy 

included  in  the  vectors  and  .  respectively.  In  case  of  Method  3,  vj, 

for  i  *  1.2....N’  will  also  be  included.  Then 

*n„™  =  S|T)  rf,:;1  "£  7..V  <6> 

where  0  f  £  1  and  7*  (Tl  ire  matrices.  Clearly  ](I)  Hi  in  the  case  if  Method  1. 


I  r.  the  cases  .4  Mith-ds  2  ar.<!  i,  has  sore  elements  equal  to  T  wnuh  muitiplv  the 

4 ^  ^ 

A  1  ^  <  o  )  components  f  21/v  (c<1  'Otherwise  0*1.  I  he  matrix 

t/#.,  (JU) 

Is  made  up  of  submatrices  consisting  of  zero  matrices  and  the  matrices  o' 


Otherwise  3 


The  matrix 
o<-  *•  * 


9 


The  linearized  boundary  conditions  on  the  k-th  iteration  have  the  form 

It  is  reasonable  to  assume  that  the  matrix  ft  t  C.b  is  nonsingular.  For  later 

*  s- 

reference  it  should  be  observed  that  the  elements  of  ,  which  multiply  the 

its  of  in  the  vector  il)  in  the  above  equation  are  all 

JRt\  (  o*t  t* 

Therefore  Ct  Of  A  t>  ( ^ 

<  I  r  f  i  n  /  R  i  i 

Clearly 


components 
zeros . 


From  equations  (6)  it  follows  that 


Since  CBQ±p  to/  -  C.  -lA  (C) 
fi  />*,  *  '/ftt  / 


I  ■>- 


}  v*  V/-L,  . 


From  ecuation  (6)  ir  follows  that 


hi 


Let  '  signify  the  Euclidean  norm  of  a  matrix.  Then,  using  the  properties 
of  norms, 

1  iT 

i<y*  «>/i  il  /s  </  -/  *) j d^/is^f  -f')U% 


3v  the  mean  viije  theorem  for  integrals. 


f 


and 


0  £•  t  h  T 


Therefore, 


where  -  T  -  > 


Hence 


(7) 
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where  c  is  a  constant  independent  of  k.  The  existence  of  c  is  predicated  upon 
the  boundedness  of  the  elements  of  (B^  +  C^)  ^  ,  C.^  ,  and  S,^  within  the  region 
of  interest.  The  elements  of  ,  C,^,  and  are  first  and  second  partial 

derivatives  except  that  the  elements  of  S,  contain  factors  At  and 

( '  )*  (  ')*  * 
y  ,  (T)  -  z (T)|.  In  Section  V  it  is  shown  that,  under  certain  reason- 

jk  j*  * 

able  assumptions,  for  any  d  >  0  it  is  possible  to  take  max  At.  sufficiently 

i  '■ 

small  that 


for  i  »  1,2,. . . and  k  •  0,1 

.  ( i )  * 


j„' 

(i)* 

(Since  z,^  is  an  approximation  to 


v,'  "  the  above  result  is  not  surprising  provided  At.  is  chosen  sufficiently 

.<  1 

small.)  A'nder  the  assured  conditions,  it  is  therefore  possible  to  choose  the 


positive  nuaper  ;  ,  wh ’ ch  contains 


.< 


as  a  factor,  as  small  in  magnitude 


•rax 


as  desirec  cy  cnoosir.g  ;  to  ie  sutficientcy  smai*. 

It  has  been  assumed  oat  the  standard  quas i linearization  method  applied 
to  equations  (3)  and  (A)  will  converge  quadraticaily  to  the  true  solution  tit. 
This  means  that  there  is  a  positiv.  contanc  £  such  that 


b  IT)  -».(T)f<£f  b  -JP  (T) 

'  Bt  I  '■  1  1  “  f?  r~ 


'3) 


for  ail  values  of  X"  in  the  interval  [0,  I],  provided  p.  is  sufficiently  close 

. .  .  . 

tc  p,  .  _ncer  tne  sar.e  assumption  , 

Let  s  be  anv  number  in  the  interval  (0,  1).  It  will  be  shown  that  if  p. 

K 

is  sufficiently  ciose  to  p^  ,  then 

implying  that  the  modified  qu.is ilineari za t ion  algorithms  are  contraction  mappings. 

-A 

Moreover ,  p  *  ?-  as  <  -  ”  ■ 

In  order  to  verify  inequality  (9),  the  triangle  inequality  will  be  employed 
as  follows: 


l/v>  y>i 1  hl/i*,‘T''A'T,l 
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By  inequalities  (7)  and  (8), 

l^\T>-£roU<:  (10) 

By  means  of  inequality  (10)  it  car.  be  observed  that  if  a  modified  method  of 

quasilinearization  converges,  then  the  convergence  becomes  quadratic  in  character 

as  c  -  0  (i.e.,  as  max  it.  •*  0). 

i  i 

Provided  p.  is  sufficiently  close  co  , 

fll?*T)-/r<Vr  -  ?  I^(T1  ~Prir)l  (llj 

For  all  values  of  T  in  the  interval  t0,  1]^ 

£  c  lF/(T)-frLr)  \*  c  \p*(r\  -friT)  I 


.axing  c  i  «e  attain 


c'A.r-Afr)l  ^  i 


tor  all  values  o:  7  in  0, 


3v  inequalities  v  10 J  - 


for  ail  values  of  T  in  [0,  l!.  The  latter  inequality  implies  inequality  (9) 


which  was  to  oe  proven. 


V.  Examination  of 


i-»ti)*  -»(i)< 


In  order  for  the  analysis  of  Section  IV  to  be  valid  it  must  be  argued  that, 

* 

for  .an?  £ :  v  #*  oc si  rive  n'isser  - 

I  -.It)*  -±(a+  I  * 


^  (T’  m  "'3 


0,1,2,...  provided  the  njxir.un  interval  size  is 
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Take  N  sufficiently  large  and  max  At,  sufficiently  small  that 

i  1 


y>^\  ic  -  t)  ~3c  (r1  -  "u 


for  i  *  1,2,...,N,  where  will  be  chosen  sufficiently  small  that  inequality 
(13)  will  hold.  It  is  reasonable  to  assume  that  inequality  (14)  can  be  satisfied 

Mi)* 

because  che  boundary  conditions  imposed  noon  z  ( r  )  have  been  chosen  in  such 

-»(i)*  °  -Mi)r 

a  way  that  z  ( T  )  will  be  a  constant  or  linear  approximation  to  v  (T  )•  The 

o  '  o 

approximation  would  be  expected  to  improve  as  max  At.  is  decreased. 

i  1 

-2k*  A 

Also  assume  that  p  is  sufficiently  close  tc  0-  that 

o  y  •  T 


where  a  will  be  selected  later. 

Mathenac ical  induction  will  he  employed.  Assume  that  it  is  known  that 
inequality  (13)  is  true  for  K  »  0,1,2,. ...k  and  i  «  1,2, ...,N.  It  will  b? 

shown  that  it  is  true  for  n  »  k  +  1.  It  was  shown  in  Section  IV  that,  under  the 

(i  )*  m.!)*'  -  * 

given  assumptions,  inequality  (9)  would  hold  true  as  long  as  max  y,  -  z.  . 

_  f 

There: ore 

for  ,\  -  0,1,2 . k . 

From  inocuaii t ies  (15)  and  (16)  it  follows  that 


if)  -  /r  (  T  J  j  t  /3  ^ 


::r  F  *  0,1,2,..  ,k  +  1.  Therefore, 

I  -»i‘>  < 

(r;-^T  it)j  i  ,3  o, 

Now  in  all  of  the  modified  methods  of  quas l i 1  near i Z3 1 ion ,  f 
_»/.)*  -»  ^  * 
00  '  ^  ^ 


is  ' :  : n e  : o ra 


-V  J  ►  . 

where  u.  -  1  (c>  and 


2. 
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Hence 


(c) 


?r 


_ r.7 

J  V  L • 


t 


tU* 

(11 


where  the  partial  derivatives  are  evaluated  at  intermediate  values  of  the 
variables.  Assum'r.t  the  partial  derivatives  are  bounded  within  the  region  of 
interest,  we  obtain 


14 


L< 


With  K  «  k  +  1  and  !C  »  k  ,  the  latter  inequality  implies  the  following  inequalities 
resnective Iv : 


^hl/r 


1  |  i  j  /  \  t* 

l]T  W\~%  ('rhyk  to\L  * 

Adding  the  two  latter  inequalities,  we  obtain 

-V'U  I  I  a,;i*>  -/Jr  I  /  ,  i  ,<• 

I  1  Ml  -  y„  i  T)  -  X  (T,1  -  *  ir)  -1 1-  ( 1  u 

using  the  fact  that  J  >■  3  '  1.  Repeated  application  of  the  above  inequality  giv 


t*1-*  l^l^uJr  uni)*'  ../,-,  „  , 

>*,  "T-l  " !  1C  tv  'I  <n  ^ (  l+MUifi  f;'  ►  '  *  '  i 


-nere: ere 


— (j* 


/*„  (T5_>,  ,T)f  -  .*  *  "L  (  \  r  u.)  ■*  5 


where  5  -  2—  ^  • 

r*'  '  C 

inequality  (17)  holds  for  D  .  T  -.1  and  i  -  1,1,..., N.  In  order  to  prove 

inequality  f  1 3 j  for  K  *  k  *■  1,  the  nucoers  £  and  i  (see  tneaualitv  (15))  oust  be 

o  . 


tnosen  sufficiently  small  tnat  z 


'1  *  -la’ 9  '1  2  ;  i.e.,  such  that 


r$!L  <-11011  *n' 

r*'.  i  -3  “  v'  '  7TV.v 
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VI.  Orders  of  Accuracy  of  the  Modified  Methods 


The  question  arises  concerning  che  orders  of  accuracy  of  the  converged 
solutions  resulting  from  che  applications  of  tne  three  modified  methods  of 
quasilinearization.  An  interval  of  lengch  At  will  be  considered.  For  simplici¬ 
ty  of  notacion  it  will  be  che  first  subinterval. 

Define  the  operators 


...  ■£  a 

Tne  .ay  lor  series  solution  at  the  dlf  terer.Cia.  equations  y  -  r  '  t ,  y tnrougn 
f our: -order  terr.s  is 


-/  r  A 


r  .u  It4;  P  -  r}iyrPpfi-  P  f~f  pf]  +  J P J 

-  -  -  /  /  c 


where  the  subscript  "0"  indicates  evaluation  at  (t^,  v^).  Expanding  further, 
we  cbcain 

t'l?  - P  ?), 

J  J 

tr  ■> r  T)  /«■  rj 

+l3p'¥?pj¥PlPp/} 

r  7  f  ctf  *  r5 ct  *■  % P  *■  ' 


(lb) 
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Now  suppose  z(c)  is  some  aporoximacion  CO  y(t)  over  the  interval 

(t  ,  t  +  At].  Suppose  y(t)  is  a  solution  to  the  differential  equations 
o  o 


We  will  think  of  y(t)  as  the  converged  solution  resulting  from  one  of  the  modified 

methods  of  quasilinearization.  The  Tavlor  series  expansion  of  y*(t  +  It)  as 

o 

determined  from  the  above  differential  equations  will  be  compared  to  expansion 
(13).  Let 

]  *  p  . 

Then  in  the  same  manner  thac  expansion  (13)  was  obtained,  we  nave 


+  +PjjJ0 

+  -r  \r  i  q  +  7  P  Q  t-  U  L,  t  4  V .  4  I 

6  Wcr  Lj/J  ft  Jj jt  j>j  jj  C 

+■  P,  Ac'[  j  *•«?/  V**  T) 

-  ^jjcr  ^  V/  /r  7/ 

T  7  7  ,  J  /  !  t  /**"  /  -1  r-  5"  j 

Jrr  ~  ijt '  Ju  '  ' 


■-here  certain  terms  nave  been  emitted  because  P ^  P'j t  ~  «-* 

'.e  terms  of  the  right-hand  side  cf  aquarians  (19)  will  be  expressed  in  terms 


;r  :  and  its  ;er ivat.ves .  Thus 


/«:  rr(t:Ju  r/cP*' 
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where 


7 


£  --  X  V 
</  j  yj  ^ 

^T-lP:j  -  ^lt'}T 

J 

r  L-7  y  fyr,  p]if'jJfL  CJJ~ 1 , 

J J J  J  y y  7/  y  / 


•-nere 


-\  ff  p  (t  y 

-7  7  7  7 


777771  7^77 

/f:r-  ":rtltyJ  *  Ptf  jj  '  */c0!  ) 

/ 

":7  vr7/~7) 

>i]-WetrV  -  77'7} 


^  \  T  / 


-r  / 


: ,  -  f  7w  7  ^ 7  7 '  ]  f  * f  ^ ( '  f  -  ? ' ) 


Djjt r  7^ 1  5 + 7, n  .7 '7 /; 

%?v  7'^'^ )  ^  7  y  71 77; 
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hf  y^/ytt‘3) 


Expansion  (19)  say  now  be  written  as  follows: 


[Htj> ^lz, ptjfl  - {  Jc'{  ft  tt,p 

+  -[  yy  y 

' -$4^1  -t$j y  f- ^  /  4 Pry (y } 

+ %  y  <■*.?> t  y  $ (r'  f T9 


*\jwr)''KTf 


fit,y)(f-pb^(t,^(c,p\ 


♦  fcJc"/  <:r;-  ji-iCZfijz.vJfzCZ/jt, 


?  y 


- 1 ,  ^  iT c  y 1  ^  ^  -  r  ^  /?  i' t  ■? j  l  a*  7  r  3  /  c  y  -T ) 

j  r  j  j  j T 


*•  -  5  ?r?  1  *■? -  f 7  j  < f- y  -  y  y  y  '3  ^  y ; 

-L3y  (r;  *)  *3 1  fj'r  j;  -  y  y][  f*;y  -  ? 

p;  j  t  y[  y  ( t.p  ^Yfj  y ' + 1  y  <yj  ^  y; 


+-i(^yrly;f(ry]r  rof-yy  +  r(jfsj 


T 
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We  aus C  now  expand  expressions  such  as  f(t^,  z^)  about  y^. 


T/r0,  *]  -  RfC(  i. 


+-  y.jjyp  '’^Py  fppfiiyp 
,,piyp  ^yp') 


77  %-ij  h  -  7 

J  *  '  1  *  /*■  /* 

7  '  “  ^  j  f  ?)+  D  T  I  r  A  )  -t-  '  £  :  *■  ~T  j  !  I  ”  (  ' ) 

■•;  *•/-'  ■;/  57‘'  J  f'e,r  LTJ  '/  :vv  '  fc  1  /<  /« 


'  1  >0  ■  ’’  "■  -  1  1 i  :  1  ,J 

.  1  '  j' 


VS^'u7vf‘  ,0yy-*^> 


A .  y.  ' 


-'*?  ;*  7  )  -  .-y 1 7 - 

■“>?  'J"‘-  •  '  /i 


-*  I  J 

A  y 


V  •  a-  •>,  ^  rf -t.  A' 


7  '  r  .  -7  ’  7  _  ,7  1 


■  <•  ■  •*  <  cr 


^  rv y^rpLypi  ^r( y  py 


l;;s . r. t _  <.x?.insion  'AC!,  -e  ,bu:r. 


Jit.-iW  ‘fp-  1 ' 7  V‘  f  7' '  7  ‘7 r-  7  f 

/  '  y  A/  7;  - / 

7^7  7  ’1  ■*"  4  j  r  “  i  a  1-7  (1-7)  r  7  !  p  ^  '  7 ! 

7 7j"l-  -  r"c  -j  jJ  1  j  r;  /•' 


7 f 


vjvsifvtK'tfcf'up  f 
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"  pij  ]  V  7 

^  ^  p/fffi+fpjte 

+ -5  yj#J 


-  t?/,f  *■><  V?  V?A‘ 

+  2  P^p  +  'l<V1j  fy '  ?  1  jc  1-  020  r s; 

/  J  j 


C~  ^  \-*, 


s  _or.sioered  :  o  ;  e  equal  to  :  ( 2 1 )  because  z  _  will  :e 
'It;  •  ' See  belowj 


Substituting  into  expansion  (21)  and  rewriting  the  expansion  through  third-order 
terms,  we  obtain 


y(tc+jc)  -ZL  ^-air  Lf-  L  J 

rciiy*rp.rj^l 

*  1  rj'fl 


j 


*  *-  \  “  Wf  T“ 


Therefore 


V.  t^t)-  7^^^  *  1  Jr  ifc 


y 


/  3  f  -  ,  ^  ^  N'**,  p 

?  Jc  l_-3(  Vjy'/  *3(VnrJf  4  r 


rr 


J  J  ' 


*«■  --*1 


-o 


^r<AC") 


In  orcer  for  t.-.e  ibove  expansion  to  as?ree  with  expansion  .13)  througn  third-order 
terms,  it  would  only  be  necessary  that 


z1  -  (  P„  £ 

l  3 

1  -J  c  L 

J  -j 

luppcse  ;7  C.  ji  r.  Then  we 

les  ire 

£ 

— *i  ^ 

/ _ ,  ^  \ 

/  j  -*)  { y  t  - 

j  r  « 

( p  *- 

- 

! 

This  is  equivalent  to  3n  ( I  -  l  )  •  1 

2 

or  in  ‘  -  h 

-  1  -  0. 

Minimi 

•n 

i 

3  i  '  -  3  i  ■*  1  with  respect  to  i  ,  we 

::nd  j  -  -j. 

Therefore , 

we  lot  V 

3 

and  we  cDtair.  tre  tern  —  (j  *  ,  •  i 

**  y  y  o  o 

n  expansion  1 

221  rather  t 

■  tin  vD  ‘  t 
V  V 

Consequent ly ,  expansion  (22)  nearly"  compares  to  expansion  (18)  through 
third-order  terms. 
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Hence,  if  y(c)  is  approximated  by  means  of  p  - ^  tC~(AC  )  , 

we  obtain  y  to  nearly  third-order  accuracy  from  modified  quasilinearization. 

Suppose  we  take  ^  f  [_  £  ♦  H  (CutAtQ  as  in  Method  1.  where 
is  the  solution  to  ■),,  -  f.  •  Then 


;• 


-- jT  r AtC'u.,p  -  sat' 


I 


*e  see  tnac  2 


r  (It ) ♦  Therefore 


"-4 
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-x  ^ 


In  the  case  of  Method  3.  t  -  yQ  +  (t  -  tQ)fo  so  chat  again  zq  -  VQ  “  0  and 
z"  -  O'.  However  z'  -  t  .  Therefore,  in  the  cases  of  Methods  2  and  3, 


-a*.  Ck 


>  'h  ^  r  c 


1  -  c  + 

■y  “  c 


r  7 


-here  £  >  I  and  3=0  in  the  cases  of  Methods  2  and  3  respectively. 
S’ow  suostituting  into  expansion  (21),  we  obtain 

1  a-tz  /  ^  t-  C  £  ) 


r,.T  W  -  "L  r  ATP.  r  3JH  rr  > 
-  .-/ 1.  <-  <• 


■J 


-  -  irV  *  -  " ) r ^  i  ir/? 

t  ■  l  Lr  lO'J.  l-r  ■/  -i 


rj 

_  -x,  -»  -i 

J-'2-£P™  +-  Ar  f  rM  ^ 

*■  )  J  i  j  L 


'J 


J 


-  JH  PCtr  ’  KJfy*  -^Yrj  Py)F 

+yy»yj)c 

+■  J f“, ,  *•  '*!,(!,  f  /( fr1 

■  -j?  /  y  y  J 

—  ->  ■  -s»  i-*  ^  1  -  ^  —  ”)  r 

1  r  _  (  fV£  )  C  C  f  iD  C  -1 —  3  J  r ?  1 
*-  •  'f  1  w ' '  u’^r  -  1  J 


Aj  y  ’Y  “  j  ?) 


[  l?,fj tl  -  L  ( R,  P, *  )  /Jir  ,'P;  f;j  H '  p )], 


.1  ) 


r  j 


because 
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4  ^  C  ^  -X  , 

where  -  *  *  +  f  £ .  4 If  a  and  b  are  constant,  cnen 

Therefore  the  expansion  becomes 

^:r -«}  =  ?.  -afF.  <-  t^-V^  A 

'  r-  ^ f  5  n 

,  />  /  _  ^  , 

~  2-^4^  r  *  ru  F)  4L  rru  E  +  ft-  FQ  f  y)? 

'  j  )  -  i  *  -!  F 

*  ^ar’/f  ■*(?'.  #3  pc*  4.  p  p)f 
<:  rfr  '7  zr)  JUJ  'J 


+  ( 


3^3  Vj^t:)(?^FJ) 

lF  J  'J  4  -  J 


+  Fjlr"tH/+‘%r}W}l/+t-uisJ 


The  latter  expansion  compares  terra-by-term  with  expansion  (IS'.  This 
that  Metnocs  2  and  3  give  rise  to  fourth-order  numerical  accuracy. 


implies 
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ABSTRACT 

The  problem  is  real-time  guidance  of  an  air-to-air  missile  engaging 
an  accelerating  target  whose  velocity  vector  can  be  predicted  as  a  function 
of  time.  A  Zernelo-cype  model,  consisting  of  the  nonlinear  translational 
equations  of  motion,  is  employed.  It  is  assumed  that  the  velocity  magnitude 
of  the  missile  is  a  known  function  of  time  on  each  guidance  cycle.  The  per¬ 
formance  index  is  a  weighted  sum  of  squares  of  the  oiss-discance  and  the 
tiae-races-of -change,  u^  and  u_.of  the  flight  path  angles.  It  is  to  be 
minimized  with  respect  to  u^(t),  u]r(t),  and  the  final  time.  A  modified 
quasilinear ization  method  is  used  to  solve  the  two-point  boundary-condition 
problem  associated  with  the  necessary  conditions  of  optimality.  It  is 
applied  to  several  missile  engagements  and  is  shown  to  be  considerably 
more  accurate  over  a  single  guidance  cycle  than  a  method  based  upon  the 
linearized  translational  equations  of  motion. 
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I.  INTRODUCTION 


There  is  an  acknowledged  gap  between  the  theory  and  the  application 
of  nonlinear  optimal  control  theory  -  especially  in  the  case  of  real-time 
applications.  Since  che  necessary  conditions  of  optimality  frequently 
cake  che  form  of  a  two-point  boundary-condition  problem  ,  che  quasilineari¬ 
zation  (or  generalized  Newcon-Raphson)  method*  comes  to  mind  because  of  its 
favorable  convergence  properties.  However,  there  are  difficulties  asso¬ 
ciated  with  the  latter  method  because  of  the  need  to  solve  sequences  of 
linear  two-point  boundary-value  problems:  One  must  deal  with  the  problem 
of  storing  intermediate  solutions  or  must  numerically  integrate  a  large 
r.umoer  of  differencial  equations.  Methods  based  upun  parameterization  of 

the  control  and  use  of  nonlinear  programming  oring  forth  similar  problem. 

2 

Modified  quasilinear izacion  mechods*  have  been  developed  in  which 
the  problems  of  storage  and  numerical  integration  have  been  reduced,  espe- 
tially  tor  problems  in  which  high-order  accuracy  in  the  numerical  integra¬ 
tion  is  uncalled  for.  As  a  matter  of  fact,  in  most  reai-cime  applications 
it  is  unnecessary  and  unwise  to  spend  computer  time  integrating  with  high- 
order  accuracy  when  the  physical  constants  can  not  be  determined  with  com¬ 
parable  accuracy. 

In  tne  modified  quas i 1 inear iza t ion  method  which  is  applied  in  this 
paper,  the  solutions  to  the  differencial  equations  are  divided  inco  several 
su'oarcs.  Intermediate  approximate  solutions  are  taxen  to  be  constant  over 
each  subarc,  thereby  permitting  the  linearized  differential  equations  to 
be  integrated  in  closed  form  over  each  subarc  in  the  case  of  the  air-co- 
nr  missiles  problem.  It  ran  be  argued^  that  this  quas 1 1  inear  izat  ion 
method  will  converge  under  certain  assumptions  and  that  the  converged  solu¬ 
tion  will  have  nearly  third-order  accuracy  in  the  sense  that  It  is  equiva¬ 
lent  to  Integrating  the  differential  equations  with  an  approximate  third- 


J 
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order  formula  of  numerical  integration. 

The  development  of  microprocessor  technology  and  improvements  in 
estimation  methods  now  make  it  feasible  to  incorporate  nonlinear  real-time 
guidance  algorithms,  such  as  che  one  developed  in  this  paper,  into  the  con¬ 
trol  systems  of  air-to-air  missiles. 

The  new  guidance  method  is  expected  to  lead  to  smaller  miss-distances 
than  the  commonlv  used  proportional  navigation  method  which  is  non-optimal 
for  the  cases  of  accelerating  targets  and  nonlinear  equations  of  motion. 

It  should  be  mentioned  that  there  is  a  system^  of  finite,  nonlinear 
equations,  involving  elliptic  integrals,  whose  solution  provides  .he  opti¬ 
mum  solution  to  the  air-to-air  missile  problems  of  the  type  considered  in 
this  oaper.  However,  the  method  is  restricted  to  two  dimensions  and  con¬ 
stant  missile  velocity  magnitude.  It  also  requires  a  good  initial  approxi¬ 
mation  to  the  solution. 

Models  jf  tr.e  physical  preblea,  similar  io  that  to  be  employed  in 
inis  paper,  nave  seen  _.seu  m  the  oast;  tor  example,  one  paper  manes  use 
o:  a  linear  expansion  about  a  solution  to  a  simple  linear  control  problem 
in  order  to  obtain  an  approximate  solution  to  the  nonlinear  equations. 

The  latter  paper  makes  use  of  several  simplifying  assumptions  such  as 
constant  closing  rate. 


II.  THE  MOD  El 

The  space-fixed  coordinate  svseec  Co  be  employed  is  depicted  in 
Figure  1.  The  angles,  >  and  <t,  define  the  direction  of  the  instantan¬ 
eous  missile  velocity  vector  V 


The  equations  of  relative  motion  of  the  target  and  missile  are 


R  -  VT(t)  -  VM(t)p(Y)  (1) 

where  R(t)  is  the  line-of-s  igh t  vector  from  the  missile  to  the  target,  V^(t) 
is  the  predicted  velocity  vector  of  the  target,  V^(t)  is  the  velocity  mag¬ 
nitude  of  the  missile, 
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Here  V  (t)  is  assumed  to  be  a  known 
adjusted  from  one  guidance  cycle  to 
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wnere  *(t;  is  a  given  continuous  positive-definite  weighting  matrix  and 
x  is  a  given  positive  constant.  The  index  J  is  to  be  minimized  with  respect 
to  u(t;  and  the  final  time,  t „ . 


:::.  necessary  conditions  of  optimality 


It  is  necessary  that  the  Hamiltonian 


*  -  u 


be  a  minimum  with  respect  to  u  at  each  instant  of  time.  Therefore,  when 
u  Is  not  on  a  boundary,  «  0  so  that 
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u  — w-1  3 


The  adjoint  equations  are 

t  »  0 


-VT>  -» 

v  o'  y 

M?Y 


(2) 


where  A  and  3  are  the  adjoint  variables  and  p  is  the  partial  derivative 

U  — Jfc 

matrix  ip/  uv  . 

The  transversalitv  conditions  are 


■  kR(tr) .  d(t_;  »  ) 

R  t_)  -  0 


'  t_  tree) 


■-.ear  .v  i  1  _  j  *  -  a  no  •  3  .<3  ‘ —  . 

Now  consider  1  tine  interval  )n  wmon  j  is  on  a  ooundary.  later 

arguments  mdic'te  that ,  it  1  is  n  a  boundary,  it  is  likely  to  oe  during 

1  T  — ^  — A  * 

the  initial  portion  ; r  tne  flitnt.  Recall  tnat  H*  *  y-  u  u  :  u  must 

>  i  minimum  with  rest-'  :  re  ...  At  ar.v  given  tiae  tne  surface  t  *  H* •  u 
dust  te  concave  rpw.irls  n-  aose  -  is  csitive  ierinite.  Let  1*  »  -V  ’c. 
~-«r.  we  nave  the  tollowmg  necessurv  condition: 


1  <  u .  ■  a  ,  then  u .  »  u 

1  ain  —  i—l  tax  1  1 


i  din. 


(3) 


;  max 


tnen  u 


t 


Since  c  and  -  are  continuous,  a*  <no  thence;  u  will  be  continuous. 

At  least  in  tne  case  or  two-dimensional  problems  with  *  1  0,  one 


wou.il  ex oec 


gr.t  pstn  1  r.  g  . 


to  varv  aonoton • o a  1 . v  with  respect 


n  tide.  This  point  will  he  argued  acre  precisely  as  follows.  If  tp  is 

.  1 

free  fn.t  given; ,  then  «  ‘  r ;  •  R't„;  *  0.  Therefore,  if  V*  (t_)  is  large 

r  r  Mr 

—*  -i  -X  -k 

in  comparison  to  V^(t_)  ,  equations  (1)  imply  p('(t  )]  •  R(t  )  *  0;  i.e. 

*  p  r  r 

the  opciaai  flight  path  of  tne  missile  is  approximately  perpendicular  to 
the  1 ine-of-s  ight  vp.  *  r  at  the  final  tine.  Now  suppose  that  at  an  in- 
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scant,  u^Ct)  *  0  for  some  T  in  Che  interval  (c^.c^,).  Therefore,  6^(1)  »  0 

— 1  — ^ 

and  equations  (2)  imply  p„  [Y, (T)]  •  R(t_)  »  0.  However,  p  is  perpendicular 

ya  a  * 

to  p.  so  that  piY  (T)]  is  parallel  co  XU^.)  and  hence  nearly  perpendicular 
A 

to  PiY  CtF) j.  This  situation  can  not  be  possible  if  v ^  varies  by  less  than 

90°  over  Che  course  of  the  flight.  Consequently,  in  such  a  case  it  can  not 

be  true  chat  u  (c )  ■  0.  It  fallows  that  u  is  monoconic  and,  since  u(t  )  *  0, 
A  A  » 

u  can  not  change  sign.  Therefore,  y  is  also  monotonic  under  the  assumed 
conditions. 

In  scnmarv,  tne  necessarv  cone  it  tons  consist  of  the  differential 


!> 

i 


i 

I 


t  r.  'i  *  r.  t?  i  ou  r.  «i  r  i 


wnere  .  is  determined  by  mean;.  ,1  •■quatlens  O). 

IV.  MODIFIED  'Il.Aa  I  LINEARIZATION  METHOD 


A  modification  of  the  quas 1 1  inear i za t ion  method  has  teen  developed 
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by  the  author'.  The  raechod  will  now  be  described  briefly. 
Consider  the  differential  equations 


2. 

y 


(5) 


wich  soae  associated  boundary  conditions  applying  at  times  t  and  t_. 

O  F 

Here  y  signifies  a  vector  of  n  components.  The  time  interval  [*  ,c,.j  is 

divided  into  subintervals  [t.  ,.t.i  where  i  -  1,2, ...,N  and  t  «  t„ .  The 

i  —  i  i  r 

linearized  equations  are 


y»:(t,v)  +  :  ( t  % v  )(y-v  ) 
s  V  '  s  s 


'  h 


.Here  v  t!  is  an  approximate  solution  to  the  twe-point  boundary-conc i c lot 
T?2C  -  oolem.  The  usua.  quas  1 1  ir.ear  izat  ion  method  tails  for  tr.e  solut:; 

to  the  :?2C  prsolea  associates  with  equations  lb;.  The  solution  y  is  the 

new  i pprox ina : ion  v  to  he  usee  in  ecuations  ;b)  tor  the  next  iteration. 

.  s 

The  solution  to  the  TP3C  problem  during  each  iteration  calls  for  several 
integrations  of  equations  hi  or  tor  the  solution  to  a  large  svstem  of 
finite  iitferer.ee  equations. 

In  tr.e  nocifieu  netr.ou  to  he  aopl  tec  m  this  paper,  the  solution  . 
to  tr.e  T?3C  prro.em  is  oDtameo  at  times  t  ,  t.,...,t  .  Then  one  lets 


.  -r  4[vt.  , )  +  v  i  t .  >  i 

-3  i.  A.  —  a  i 


n  tr.e  i-th  time  interval  for  i  *  1,2 . N.  As  we  will  see,  the  resultm 

met  rod  converges  ir.  the  missile  encounters  to  he  msiderec.  Moreover,  it 
la  cnowr.  that  the  numerical  accuracy  of  the  converged  solution  Is  nearly 
of  tr.e  third  order . 

Over  earn  >  bi.iterval,  v  in  equations  'h  will  be  constant.  Only 

s 

v  •  t  )  y  f  t .  ,  . . .  ,y «  t,..!  must  he  stored  after  each  Iteration.  Moreover,  the 

integration  of  equations  (b)  over  each  interval  has  been  simplified.  In 
face,  it  can  often  be  carried  out  in  closed  form. 

Assume  Chat,  with  appropriate  substitutions,  equations  (5)  take 


an  equivalent  form 


y  *  f (y)  +  g(t) 


Then  on  each  iteration  equations  (6)  will  take  the  form 


y  »  A  .  y  +  b  .  (t) 
1  i 


(7) 


on  the  i-th  subarc,  where  A,  is  a  constant  matrix.  It  can  be  shown  that 


y ( t  )  »  3 .  y ( t  .  )  +  i. 

1  A.  1  **  1  1 


(8) 


or  some  matrix  3  and  some  vector  Let  »  I  and  M.  »  M.  ,  3. 

or  i  1  S-l.N-l . D.  It  can  easily  be  shown  chat 

y  <.  t  )  *  M  yi  t  )  *  M.J.  . 


Therefore,  using  the  latter  expression  which  gives  v(t_)  in  terms  of  v(t  ), 

'  o 

'ne  :ar.  write  ail  the  boundary  conditions  in  terms  of  v(t  ),  One  deterTiir.es 

o 

y  t  ;r.  : r ; e r  to  satisfy  the  boundary  conditions  arc  eoploys  equations 
•  m  :rder  to  obtain  y  <  t ,  _■  ,  y < t y ( t„) . 

here  we  have  assumed  tr  is  given.  However,  it  is  possible  to  adjust 
it  in  an  "outside  loop"  as  will  se  explained  later. 

The  quasi  1 ..near i gat  ion  method  can  also  be  modified  <s  follows:  Rather 
than  accenting  the  solution  y  to  the  linear  tied  equations  as  the  new  approximate 
solution  y  ,  let 

-*  'new;  -»  'old)  .  -»  fold! 

1  t  .  ;  »  V  '  t  .  )  ♦  l,v(t  -v  ■  t  , 

S  1  '  S  1  IS  i 


!  .r  .  «  ,,N,  where  i  is  a  posit!’,  e  constant  chosen  such  that  the  maxi- 

m  .m  n.ihgc  in  y  will  not  he.  me  "dangerously"  large.  This  is  analogous  to 
the  damping  method  sometimes  used  in  connection  with  the  Newton-Raphsor 
method  used  tn  solve  a  system  of  nonlinear,  finite  equations.  The  damping 
technique  has  enabled  the  qua sil lnesr i za 1 1  or  method  to  converge  in  some 


cases  in  the  air-to-air  missile  problem  when  it  would  not  have  converged 
otherwise.  (In  some  cases  it  might  be  desirable  to  take  a  >  1,  but  we  have 
not  done  so.) 

Now  we  will  apply  the  modified  quasilinearization  method  to  the  air- 
to-air  missile  problem.  Although  it  is  applicable  to  che  three-dimensional 
problem,  it  will  be  applied  herein  only  to  the  two-dimensional  problem  with 

1  t  .-»  2  1  ,:f  2 

v„  1  0.  let  v  »  v  Take  J  *  —  k,R(t_)'  +  —  /  -a  dt,  where  u  *  y .  Eaua- 

-  A  r  t.  o 


i 


>9) 


-1ST- 

wnere  p  ‘  •  (cos  -.sin  v)  and  -  in  orcer  to  place  the  ecua cions  in  che  form 
it  equations  ,7;  -  w:  introduce  tne  cilferencial  equation  V  w  qtt)  with 
q • t .  being  a  given  function. 

The  linearized  equations  on  the  i-th  subarc  are 


where  all  quantities  with  the  subscript  i  are  constant  over  the  i-th  sub- 

-»  a 

interval.  The  quantity  p  ,  for  example. stands  for  p  (y  t.  Thus  the 

y  si 

subscript  i  signifies  evaluation  on  the  approximate  solution  v  on  the 

si 


i-th  interval. 
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Lecting  r .  »-V  p  a  -  V  p.  A.,  and  taking  advantage  of  the 
i  Mi  Yi  l  mi  1 

knowledge  that  A  is  a  constant  and  p  »  -p,  the  system  (10)  can  be  written 


To  i  o 


1  M  ! 

i  _kT_k 


;  u  ;  »  |Oi  0  o  u  j*  ;V+[VMi'v«(t)1^iV',i,i  ! 


o  oi  s 


■VX(t)pi4VT(t,’'i'i 


.  .e.,  on  tne  i-th  suointervai 


y  *  A .  y  +  b  .  ( t ) 


Let  .  Let  cs;  =*  cos(-';lt),  sn  .  *  sin(-'JLt),  and  i.  ■ 

in  the  case  of  i  _  0.  Let  :s  =  cosh  ( ■' .  it ; ,  sa.  «s  inn  .  it  \  ana  i. 

-hen  a,  ■  j.  Here  we  are  assuming  a  constant  interval  size  it.  Let 

.  3 


•  ■  .  *  3  D  .  1  ;  .  ■  ,  .  * 
>  -  l~:s  .  -3 :  t, .  Lt  sn 


*  •  .  _t-sn  .  ;  %  . ,  ana  .  .  » 

1  1  1  -*  i. 

When  \  it  is  snail,  it  is  necessarv  to 


evaluate  .  ,  - ,  ,  . ,  ano 


3y  -leans  or  series.  For  examne,  vnen 


'  i -  1  ,  is  *  1  ir.ii  -  m  ;  r  le  r  to  ivoid  suo traction  or  nearly  ecua*  nuir.ce 

one  expar.es  is  .  me  evaluates  -  ,  .  using  tne  in:  mite  series 


it  :  . 


*“  k*0 


-  ,  .  .  tw  ,  -  .  ^  >  1  k. 

rrorr  -v  '’i-° 


(2  r 


in  tr.ts  formulation  it  is  assumed  that,  on  eacr.  subinterval  VM  f!  * 
Vt,_,  '  *  ;  '  t  -  c  ;  _ .  )  /  ^  t ; :  V  .  t . ;  -  V  '  t  .  _  .  )  ;  and  similarly  for 

It  :ar,  be  shown  basing  Laplace  transforms  or  matrix  eigenvalues  ana 
eigenvectors;  treat 


R(t . ; 


>  !  t  .  I 

1  “  w 


■iU  ,  ) 


f  - 


* 

0 


vhe  r  e 
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I  csi 


B* 

l 


'u!  .  , 

i,  C.  -  IjIVjj.SjVjj]  , 


.  5 .  '  .  an ,  cs . i 
ill  lj 


Vli  I 


'31  1 


Pi  .3.v  -o . ; .  , 

i  i  2i  i  i  3i 


,  0.  -  |5.c  v-.-C.:;.-* 

i  I  i  i  2i  i  i  11  i 

l  ^  j 

.:(;lv4i-5iYlV3i)^+&t(ti*-cc;l*/2) 


— — 'V  ( t  1 -V  ( I  , )’p  "  \  , 

If  M  1  M  1-1  '^1  l 

1  _>  i 

■  -  -  y  m  -  ~  '  v  v  •  ~  :  _  •  •'  1  r  : 


2  ’i'Vi  1 


1  1 


«'e  -ill  1c;  the  M  matrix  r:  the  general  formulation  be  ’-Titter,  in 
terms  of  i:<2  ratrites  as 


.et  M,.  3  1  anc  Nv  -  0.  Therefore 


*>  *  » 

M  .3  ,  N 


#  *  # 

SB  -el 


t,S-i , . . . , a.  Then  the  approximate  sc-atlan  to  the  linearized 


enua  1 1  :ns  1  is 


■  r : 


r 

u  (  t 
•"  <  t_; 


u  ( 0)  + 

_ * 

R 


i-1  N 


-*  *  -a  t  -u,  --  * 

(„  ^  V  -t-  <b  )  . 

r  1  1  1  ‘ 


In  order  to  solve  the  tvo-point  boundary-value  problem,  one  must  set 


u(tj.)  »  0, 


kiiftp),  and  must  solve  the  above  svstem  (of  linear  equations) 
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for  u(0)  and  R(tf) .  Then  one  calculates  the  values  of  y ,  u,  and  R  at  other 
times  using  equations  (11). 

Here  it  has  been  assumed  that  t^  is  given.  However,  in  practice  we 
have  starred  with  an  initial  approximation  to  t^  and,  after  the  iteration 
with  quasiiinearization  seems  to  have  stabilized  somewhat,  we  have  corrected 

— k.  — X 

t_  after  each  iteration  using  the  condition  g(t„)  “  0,  where  g  »  R(t  )  •  R(t_). 
r  r  F  r 

We  let 


(New) 


(Old) 


with  g  being  computed  only  approximately. 

No  bounds  have  been  placed  upon  u  in  the  numerical  examples  to  be 
presented.  However,  the  procecure  can  be  extended  in  order  to  induce  me 
:ase  of  boundec  u.  The  sub  intervals  may  be  ad'usted  fusing  the  function 
u*  of  conditions  .3;  in  a  manner  similar  to  that  used  to  predict  t_  from  g) 
in  order  that  a  will  not  be  rounded  over  only  part  of  an  interval.  On  each 
interval  on  wnicn  is  on  a  oour.carv,  the  second  of  equations  (9)  would  be 
r<  p.aced  r.  i  «  1.  Equations  11  •  ,  etc.  would  have  to  be  modified  accors- 


V.  A  LINEAR  METHOD 


for  the  same  of  comparison,  we  will  develop  a  rather  elaborate  linear 

method  for  optimizing  the  air-to-air  missile  problem. 

Let  v  1 1;  ‘  an  initial  approximation  to  u  solution,  where  v  (t)  ■  v 

S  S  SI 

?  constant  on  tne  i-th  sub lr.tervc '  .  The  linearization  of  equations  (1)  on 
the  i-th  interval  is 


) 


R  - 


vT(t)  -  vMft)p3  -  vM(t)Py; 


where  l  -  1,2,. 
is  to  choose  u 


iy  which  minimizes 


P(. '< 


and  p 


ys 


P.. 


). 


The  idea 
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1  ,2  1  y  t  -»-» 

J  *  f  klR(t_)r  +  j  f  ( Y  +U)  W(t)(-r  -hi)  dt  . 

Z  f  z  o  s  s 

A  development  similar  to  the  derivation  of  the  necessary  conditions 
of  optimality  of  the  nonlinear  problem  gives  the  necessary  conditions 


8  *  kVMpys  R(CF} 

^  ,,-l  - 

y  »  -*  o 


llearlv 


s(tp)  -  0 


d(t)  »  k(.\C  dt)R(tr) 

"F 


V  t )  »-kW  *(/'  V  p  ‘  atlR(t) 
t_  M  ys  r 


v  u  t 


Y(t)  -  v  -  kC(t,  t_)  R  (.  t_) 
o 


Theretore 


'm?.s 


r  -k  -j>  *r_k  t  j. 

Rft„)  ■  R  +  ;  fV  -V  p  V  p  "(  )dt  V  p  v  dt 
o  o  T  M  s  M  >s  s  t  M  vs 

o 


Rftj.)  *  5  -  k[/,  VM?-,S  C(t,  tp)dt;  Rit?) 


where 


s  •  R  +  / /  [V,  -  7  p  *  V  p  ‘  (y  -y  )]  dt 
o  t  »  Ms  Mrs  so 


so  chat 


/f 


U  *  k  /tr  V  p  C(t.t  )dt]  R ( t F )  -  S  . 
o  .  r 

The  latter  linear  equations  must  be  solved  for  R(t^.)  in  order  to  obtain  the 
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desired  solution.  This  is  a  three-dimensional  formulation  and  is  based  on 
the  assumption  that  tf  is  given.  Observe  that  V^,  V^,  and  W  may  be  given 
functions  of  time. 

VI.  NUMERICAL  RESULTS 

This  section  contains  numerical  results  for  the  two-dimensional  en¬ 
counters  now  to  be  described.  At  the  initial  time,  t  ”0,  the  missile 

o 

is  located  ac  the  origin  and  the  target  is  914.402  meters  (3000  ft)  away 
on  the  positive  x-axis.  The  missile  velocity  magnitude  is 

V  *  217.712  t  +  295.549  m/sec. 

M 

There  are  no  bounds  upon  v.  The  target  makes  a  9g  turn.  Its  velocity  mag¬ 
nitude  is  V„.  »  295.352  a/' sec.  Specifically, 

— *  0  .  ^  0 

V_,  -  V,,  cosutj,  V  *  V,  sin(.t) 

si  L  i  -  i 

wnere  j  *  17 . 121 7° sec  .  We  will  consider  four  cases,  corresponding  to 

OOO  0 

v  *  -30  ,3  ,45  , 30  .  3v  trial  it  was  found  that  the  weighting  factor 
x  -  .00012915b  raa*-sec,n“  gives  reasonable  results  in  all  four  cases. 

A  starting  value  of  tr  -  2  seconds  was  used  in  all  cases. 

No  full  simulation  of  the  guidance  methods  has  been  carried  out. 

It  is  realized  that  a  final  choice  of  method  for  a  particular  missile 
must  depend  upon  the  results  of  a  complete  simulation  of  a  wide  range  of 
encounters . 

The  initial  ^starting)  solution  for  both  the  linear  and  quasilint  a1- 
guidance  methods  makes  use  of  the  knowledge  that  y(tp)  ”  0  and  that  y 
is  monotone.  It  is  assumed  that  V  2  constant  and  that  y(t^)  is  approxi¬ 
mately  equal  to  the  angle  between  the  positive  x-axis  and  a  line  connect¬ 
ing  the  origin  to  the  position  of  the  target  at  time  tp.  This  angle  was 
taken  to  be  40°  in  all  cases.  The  constant  y  is  chosen  such  that  y(tp)  •  40°. 
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In  mosc  cases  this  initial  solution  gives  a  value  for  Y(to>  which  is  rather 
close  tu  that  obtained  from  the  converged  quasilinearization  method.  In 
all  four  cases  the  lacter  mechod  converged  with  N  *  15  when  it  was  started 
with  the  solution  just  described. 

Figures  2  and  3  depict  the  initial  solution  and  the  converged  quasi- 

iiaear  for  the  tct' s  of  "  ■  -30°  ->nd  '■  «  90°.  Tho  order  two  cases 

o  o 

are  similar.  A  value  of  N  -  15  was  employed.  The  final  time  t^  as  well 
as  v(t)  converged. 

Table  1  shows  the  effecc  of  the  number  N  of  subintervals  for  the  case 

of  y  ■  0°  as  an  examole.  There  is  little  change  in  v(0)  and  R(t  )  as  N 

o'  r 

increases.  However,  fewer  iterations  are  required  as  N  increases.  It  has 
oeer.  found  in  ocher  cases  that  there  is  a  high  risk  of  nonconvergence  when 
N  <  15.  The  iterations  were  carried  out  until  R. (t_)  and  R,(t_)  seemed 
to  have  converged  within  about  one  meter.  Calculations,  however,  were 
carried  out  in  terns  of  feet,  radians,  and  seconds. 


Table  1.  Quasilinearization  results  for  v  -  0 


N 

6 

10 

15 

20 

v(0)  (dcg/sec) 

11.6 

42.9 

43.6 

43.6 

R ,  ( t  „)  { me  t  ers  ; 

-2.1 

-2.5 

■2.6 

-j.Q 

(Tnctersj 

5.5 

5 .  } 

5.3 

5 .  b 

Nc  .  Iterations 

19 

15 

11 

u 

16 


Table  2.  y(0)  (deg/sec), 

comparison 

of  linear 

and  quasilinear  methods 

YQ  (deg) 

-30° 

0° 

45° 

90° 

Linear  guidance 

54.5 

36.3 

-3.7 

-48.2 

Initial  solution 

70.0 

40.0 

-5.0 

-50.0 

i  f  or-3  t  i  on ) 

90.7 

45.3 

-10.1 

-68.3 

Quasi.  (10th  iteration) 

7S.2 

43.7 

-11.2 

-56.1 

Quasi,  (converged) 

72.1 

43.6 

-11.2 

-51.3 

In  Table  I  the  lar.  _c  and  quasilinear  methods  with  S  -  15  are  compared 
for  all  four  oases.  The  linear  guidance  method  is  started  with  the  same 
solution  as  the  quasilinearization  method.  .  .  seuid  be  observed  that  the 
value  ’(t  )  of  the  guidance  conroana  obtained  from  the  linear  method  differ 
appreciably  from  that  of  'he  quasil inearizat ion  method  which  gives  the  optimal 
solution  to  the  nonlinear  translational  equations  ci  motion. 

It  was  found  that  when  k  was  increased  to  .3001 5.  tne  quasilinear  itatior. 
algorithm  did  not  usually  converge.  It  can  be  seen  that  the  value  . C00012S 106 
used  for  k  led  to  miss-distances  which  may  be  a  bit  large  for  some  purposes. 

If  one  wishes  to  decrease  the  miss-distance  it  is  possible  to  increase  k 
gradually  during  the  iteration. 

A  constant  subinterval  site  It  was  used  during  each  iteration.  However, 
since  v  is  much  larger  at  the  beginning  of  the  flight  than  at  the  end,  the 
numoer  of  intervals  required  can  be  reduced  if  the  interval  sice  is  varubie. 
As  a  rule  of  thumb,  the  subinterval  should  be  sufficiently  small  that  y  does 
not  vary  by  more  than  10°  over  the  interval. 

An  alternate  method  for  .omputing  t ^  would  be  to  solve  tor  it  simuitan- 

,  -A 

eously  with  y( t^)  and  R^.,  rather  than  treating  t^  separately.  This  alternate 


17 


way  has  not  been  studied. 
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